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Strong coupling series for QCD at finite temperature and density 
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We discuss the use of strong coupling expansions for Yang-Mills theory and QCD at finite 
temperature and density. In particular we consider the onset of temperature effects for the 
free energy and screening masses, derive the hadron resonance gas model from first principles 
and compute the weakening of the deconfinement transition with chemical potential. 

§1. Introduction 



Strong coupling expansions were among the first tools to study lattice gauge 

r^ , theory and produced analytical insight into confinement and glueball masses. Here 

we discuss some recent extensions of these techniques to finite temperature and 

density.'^ Strong coupling expansions are complementary to those in weak coupling 

a , and yield convergent series with a finite radius of convergence. They are the only 

^ analytic tool to address the confined phase of QCD from first principles, allowing 

fl^' to study the onset of finite temperature effects at low temperatures. Furthermore, 

,^ , with strong coupling series of sufficient length we may hope to establish a connection 

between QCD and models resulting in the strong coupling limit, which are often used 

for studies of nuclear matter at finite baryon density .'^''21' 

> 

^f^ ' §2. Equation of state and screening masses for Yang-Mills 

''^ \ Starting point is the YM partition function using the Wilson action 



Z= I DU exp -Sym [U] , Sym = E ^ ('^^ ^ + ^^ ^^ " 2A^) > (2-1) 

with the lattice gauge coupling j3 = ^. An expansion of / = — -i In Z in /3 proceeds 
/\ • by expanding in group characters Xr{U) followed by a cluster expansion,^ 

/ = -61n co(/3) - ^ E ci{C)^cp{x,r, ^{X,) = f DU J] dr,ar^Xr,{U), 
c={x"') i peXi 

(2-2) 
where H = V -Nt is the lattice volume, dr and a,, dimension and expansion coefficient 
of representation r, and cq is the expansion coefficient of the trivial representation. 
The standard expansion parameter to express results in is conventually the coefficient 
of the fundamental representation, u = af. The combinatorial factor a{C) equals 
1 for clusters C which consist of only one so-called polymer Xi. The (vacuum) 

typeset using VTpT^.ds (Ver.0.9) 
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Fig. 1. Left: LO graph for the free energy. Right: LO graph for the screening mass. 



quantity in Eq. ()2-2p is customarily called a free energy density, because the path 
integral corresponds to a partition function if one formally identifies /3 with 1/T. Here 
we are interested in a physical temperature T = l/(aNt), realised by compactifying 
the temporal extension of the lattice. The physical free energy is then obtained by 
subtracting the divergent vacuum free energy, 

f{Nt,u) = fiNt,u)-f{oo,u), (2-3) 

and the pressure is P = —f. Group integrals are evaluated using the formulae 

dUXr{UV)Xr{WU'') = XriVW), J dUXr{U) = Srfi- (24) 

Due to the latter the contributing graphs Xi have to be objects with a closed surface. 
The graph contributing to the leading order of the free energy density is a tube 
of length Nt with a cross-section of one single plaquette. Fig. [TJ Summing over all 
such graphs on the lattice, their contribution is 






(2-5) 



Thus, the strong coupling limit at /? = (and thus T = 0) has zero free enery density 
or pressure, as expected. Moreover, we see that the pressure rises very slowly with 
T as the leading order starts at a high power in the expansion parameter only. 
Corrections to the free energy density through order u^ with respect to the leading 
order term for various Nt have been calculated for 811(2)21 and for SU(3).'^ 

Next, let us turn to screening masses. The one with the quantum numbers of 
the lowest lying glueball is defined by the spatial correlation function of plaquettes, 

C{z) = {TrUpM TrUp,{z)) = N^^J-^^ln Z{P,Pih) 



d^xdh 



(2-6) 



/3i,2=/3 



At zero temperature the exponential decay is the same as for correlations in the time 
direction, and thus determined by the glueball masses. At finite temperature, the 
LO graph for the difference to the vacuum mass is shown in Fig. [1] and gives 

Am^^[Nt, u) = m{Nt) - m{oo) = --Ntu'^^''^ (2-7) 

The finite T effect is to lower the screening mass in the confined phase and is sup- 
pressed by high orders in the strong coupling. This explains why in the confined 
phase it is close to the vacuum glueball mass as observed by Monte Carlo.'^I' 
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§3. QCD with heavy quarks 

Wilson fermions can be included by means of an expansion in the hopping pa- 
rameter Kf = {2am f + 8)~^. Defining the usual hopping matrix Mp* Grassmann 
integration over fermion fields gives to leading order in Hf per flavour 

- ^.^ = E T *'^'^^) = "^^''^)'^' ^ (e^^*L(x) + e-^^*Lt(x)) + . . . (3-1) 

I X 

Kronecker deltas in the hopping matrix M force the sum to extend solely over closed 
loops on the lattice. Again the finite temperature effects are in the difference between 
finite and infinite Nt, and hence in the loops winding through the temporal boundary. 
Thus the leading order heavy fermionic contributions to the effective action can be 
written as in the last equation, which only holds for finite T, and the dots represent 
loops that wind more than once. 

Now we calculate again the free energy density, this time performing character 
expansions in both U and L. Expanding all terms up to 0(k^ *) and doing the 
group integrals we get for two flavours u, d 

P = ]v^{4(2^«)''^' + 8(2k„2/.,)^^ + 4(2Krf)2^*| 

+ 6[2k,(2k^)2]^* +4(2Kd)3^4(e3»^^* +e-3«^^*) (3-2) 

If we now recognise the hadron masses to leading order hopping expansion as 

Mesons: arrirp = —\n2Kf — Inlufi (3-3) 

Baryons: amfff = —\n2Kf — ln2Kj/ — ln2Kf", (3'4) 

we are able to rewrite this as 



*" I 0- 1- J 



P = 

Nta 

+ J^.[^E--'^^'^^"'+^E--'^^'^^"'\^osl.{,sN,), (3-5) 

2 2 

which is nothing but the pressure of a hadron resonance gas. We have thus derived 
from first principles that the latter arises as an effective theory for QCD in the strong 
coupling regime. Note that this is a generic feature holding also for Yang-Mills theory, 
which can be represented as a glueball gas.l^'l^ 

§4. The deconfinement phase transition 

It is well known that the deconfinement transition is first order in the pure gauge 
limit, K = 0, and weakens for finite quark masses until it disappears at some critical 
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Fig. 2. The critical quark mass mc{fi^) for TVt = 1. Left: n = 0. Right: Nf = 3. Error bars 
represent the scatter for different Pade-approximants. 



Kc, which represents an upper bound on the radius of convergence for our strong 
couphng series. One may extract this quantity from the Polyakov loop susceptibility, 



xl{J) 



1^ 
VdJ^ 



Z{J)\ 



j=o, 



S{J) = SYM + Sg + jY,iL^ + Ll) (4-1) 



Again we have to perform a double character expansion. The leading graphs are 
neighbouring Polyakov loops tiled with plaquettes and their decorations. This leads 
to a double series for xl{u, ^)- At the critical parameter values Uc,hc, the suscep- 
tibility will diverge with a critical exponent 



XLit) 



1 



{tc-ty 



(4-2) 



with the parametrisation u = n ■ t,h = - ■ t. We now model the series for the Dlog 
of XL by Fade approximants in t, which have simple poles at the critical parameter 
values with the exponent A as residue. We find A = 1.03, which signals 3d Ising 
universality with 7 = 1.237. Fixing to the exact value of the exponent, the corre- 
sponding critical parameters and masses can be extracted more accurately. Note 
that Kc ~ O(10~^) and hence leading order in the hopping parameter is good. The 
results for the critical masses for a coarse A'^t = 1 lattice are shown in Fig. [2j In 
particular it is possible to compute the change of the critical mass with real and 
imaginary chemical potential. This establishes that in QCD with heavy quarks the 
deconfinement transition weakens with real chemical potential, in full accord with 
numerical findings from the Potts model.'^ 
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